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& ■ !• Introduction 



In this work, we are interested in the homogenization of energy quantities for elec- 
tromagnetic waves in the high frequency limit, and more particularly for Maxwell's 
equations. Our interest is also in dealing with interactions with various boundary 
jjj ■ conditions. For this purpose, we use a scaled variant of H-measures (see L.Tartar 

or P.Gerard), known as semi classical measures or Wigner measures, introduced in 
[14], [19], [20]. 

One of the most important predictions of Maxwell's equations is the existence of 
electromagnetic waves which can transport energy. 

For this purpose, the Theory of Radiative Transfer was originally developed to 
describe how light energy propagates throught a turbulent atmosphere. This the- 
ory can applied to representative problems involving reflection, transmission, and 
diffraction in both homogeneous and inhomogeneous media. 

Justification of this theory in high frequency limit, as well as for other waves equa- 
tions, can be found for a deterministic medium in the works of P.Gerard [10], [12], 
and C.Kammerer [15] as well as by P.L. Lions and T.Paul [18] and L.Miller [16] and 
G.Papanicolaou [21]. 

Our purpose in this paper is to describe this energy propagation for Maxwell's 
system, coupled with various boundary conditions, and with a typical scale which 
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is played here by the frequency. 

We shall consider Maxwell's system, with electric permeability e, conductivity a 
and magnetic susceptibility rj, in the half space (x 3 > 0) of R 3 , with the courant 
variable X — X ^ X ^ X 3 . These quantities arc 3x3 matrix valued functions of x. This 
system is given by the following equations 



(1.1) 



i) d t D £ (x,t) + J e {x,t) = rotH £ (x,t) + F £ 

ii) d t B £ (x,t) = -rotE £ (x,t) +G £ {x,t) , 
Hi) divB £ (x,t) = 

iv) divD £ (x,t) = g £ (x,t) 

where t e (0, T), and E £ ,H £ ,D £ , J £ and B 6 are the electric, magnetic, 
induced electric, current density and induced magnetic fields, respectively. Morever, 
p £ is the charge density (a function uniformly bounded in i 2 (R 3 ), and where F £ , 



G e a 2 (l 3 ) 3 are given. 

We complete this system by the following constitutive relations 
f 1) D £ {x,t) = e(x)E £ (x,t) , 

2) J £ (x,t) = a(x)E £ (x,t) , 

3) B £ {x,t) = i){x)H £ {x,t) . 



(1.2) 



We shall only be interested in time harmonic solutions of this system and in the 
high frequency limit. For that purpose, we look for solutions in the form 



' D £ {x,t) 


= D £ (x)^{expif} , 


H £ (x,t) 


= H £ (x)^{expif} , 


« J £ {x,t) 


= J £ (x)3?{exp^} , 


B £ (x,t) 


= B £ (x)n{expif} , 


E £ (x,t) 


= E £ (x)n{cxp I -f} , 



(1.3) 



where cj is the given fixed frequence, that we assume different from 0. Note that we 
use the same letters on both sides of the above equations to simplify notations. 

In this work, we assume that the matrix e , r? ,a, are 3x3 "scalar" matrix valued 
functions given by 

(e{x) \ 



e = e(Id) 3x3 = 



e(x) 
e(x) 



(1.4) 
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and 



fj = r?(Id) 3x 3 



/V(x) 
rj(x) 


V 



\ 





<r(Id) 



3x3 



(a(x) 

a{x) 

cr(x) 
V 



(1.5) 



where e ,rj ,a, are smooth (scalar) functions in C 1 (K 3 ). This usual assumption could 
be certainly relaxed, but at the expense of much more complex spectral calculus. 
With the above notations, the time harmonic form of Maxwell equations are then 

rotE £ - iujr)H £ = F £ , 

(1.6) 

G £ . 



IOtH £ 



<tE £ 



Note that we have not written the third and fourth equations appearing in system 
(1.1), since in fact we assume that the right hand sides of (1.6) do satisfy the usual 
compatibility conditions. 
Set 



E £ 
H £ 



El 
El 
H £ 

HI 



and 

*'-(**V- U Q o>-U Q o') ^ 

where the constant antisymmetric matrices Qk , 1 < k < 3 are given by 

/0 \ / 1\ /0 -1 0\ 

Qi = -1 , Q 2 = , Q 3 = 1 . (1.9) 
\0 1 / \-l 00/ \0 0/ 



Above, the matrix C is given by 



while the right hand side is f £ = [ qe )■ 
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Assumed uniform boundedness and symmetry of the permeability and susceptibility 
tensors show that system (1.6) is a symmetric hyperbolic system as follows 

T*+t*%: + c*-r. (1.U) 

j=l 3 

As a first boundary value problem, we shall consider system (1.6) or equivalently 
system (1.11), posed in a domain, that we choose to be the upper half plane R+ = 
{x,x 3 > 0}, with a perfect conductor boundary condition at the flat boundary 
x 3 = 0, i.e 



n+ AE £ = {) on x 3 = (1.12) 







where ii = | J is the unit outward normal vector to R 3 . Xol <( ha; i is < L< >i : iai n 

is not bounded, but this is unimportant since we will localize all our functions. 
The second problem dealt with in this paper will be a transmission problem. To 
simplify the exposition, we will consider a medium, made of two parts: R+ = {x € 
R 3 , x 3 > 0} will be the exterior medium, while R_ = {x E R 3 , x 3 < 0} will be the 
interior one, each caracterized by distinct electromagnetic coefficients. 
We suppose that our electromagnetic field is created by an incident wave u mc = 

jjinc 

In R 3 , we consider the following exterior problem, characterized by the dielectric 
coefficients (e ext (x) ,n ext (x)) belonging to C 1 (R 3 ), and scalar valued, see [17] 

TQtE ext, e _ i(Jr] ext H ext, E = q 
TQtH ext,e + iuJe ext E ext,e = Q ^ 
/ e extJ7jext, S _ '^ext Jjext, S ^ n + 



^ C 



Here E ext ' £ } H ext ' e are the so called exterior fields, r — |x|, x = (xi,X2,x 3 ) and 
the third equation is the classical Silver Muller radiation condition, see for more 
details [8], [17], with n + being the unit outward normal vector to R 3 . 

In R 3 , we consider the following interior problem, which is characterized by the 
dielectric coefficients (e mt (x) ,ri mt (x)) belonging to C 1 (R 3 ), and scalar valued, see 
[17] 

IQtE int,e _ i(J7] int H int,e = Q ^ 



rotH mt - 6 + iue mt E mt < £ = 0. 
We impose the following boundary conditions (Caldcron condition) 

E l £ nt An- - {Ef l + E inc ) An- =0 on x 3 = , 

H l £ nt An- - (H§ xt + H mc ) An- = on x 3 = 



(1.14) 



(1.15) 
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where n~~ is the unit outward normal vector to M_ . 

£int, e ^jjmt,e are tne g0 ca n ec [ interior fields. Note that there is no condition 
at infinity in the interior problem, mainly because we have assumed intuitively a 
localization near = 0. 

In the third and final part, we generalise these two cases, and we study the curved 
interface case, where the plane xz = is now replaced by a curved interface, in the 
spirit of the work of Gerard and Leichtman [13]. 

More precisely, we consider Maxwell's system (1.6) given above the surface given 
by T : x 3 — 4>{x ), where x — (x±, x 2 ), and 4> € VF 2 (]R 2 , R) is a scalar function. 
We consider this system in time harmonic form, in the high frequency limit, and 
we consider a perfect boundary condition on T. 

For each of the above cases, we shall study propagation of energy like quantities, 
using the framework of semi classical measures. Basic facts about these tools are 
recalled in Section II, refering the reader for more details to [11], [19]. 
Then in Section III, we consider the above cases of Maxwell's equations, with diffcr- 
ents boundary conditions, and in particuliar, we prove therein the following results 

Theorem 1.1. Perfect conductor case Consider time harmonic Maxwell's sys- 
tem in the half space x% > with a perfect boundary condition, written in the form 
(1.11), with solution vector u £ . Let 9{x) be a test function with compact support 
that is equal to one on a compact set K C M 3 . Let u 6 ' £ = 9u £ be uniformly bounded 
in I/ 2 (M 3 ), with (up to a subsequnce) an associated semiclassical measure jl. Then 

the semi classical measure jl is supported on the set (x € Supp 6,k € M 3 ) 

U = {(x, k) , lu + = lu} U {(a;, &),<*>_= lo} (1.16) 

where v(x) = = is the propagation speed. Above ui n = Uq{x, k) = 0, uj + = 

y/e(x)n(x) 

uj + (x,k) = v(x)\k\, ui- = U-(x,k) = — v(x)\k\ are the eigenvalues (of constant 

3 

multiplicity two) of the dispersion matrix L(x, k) = '^^(A°)~ 1 kjA J ' . 

The semi classical measure jl(x, k) has the form 

{fl(x, k) = u\(x, k)b]_(x, k) <g> b^(x, k) + a\{x, k)b\{x, k) <g> 6+ (x, k) 
(1.17) 
(x, k)bl (x, k) <g> b 1 * (x, k) + pt 2 _ (x, k)b 2 _ (x, k) <g> b 2 *{x, k) 

where n\ , M+ are ^ wo scalar positive measures supported on the set 
{(x, k) , uj + = oj} and a\_ are two scalar positive measures supported on the 

set {(x, k) , uj- =ui}. b\_ (resp. b 1 ^ ,b_) are two (normalized) eigenvectors of 
the matrix L(x,k), corresponding to the eigenvalue lo + (resp. u)-). 
Furthermore, the scalar measure p} + satisfies the following transport equation 

VfeW+.V^V - V x UJ+.V k ll\ = Ufe3fc3[^a+ 7 l*fe 3 =fc- + ^+ T l (5 fe 3 =fe+]^3=0 (1-18) 
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where v\ + , are scalar positive measures associated with the semiclassical mea- 
sure v corresponding to the boundary term u £ >®(x',Q), with x' — (£1,2:2) and 
k = k/ | k \. The wave vector A: ± (fc') = (k',kf) is defined by 



Finally, we have denoted by T\ the operator defined as follows: for all smooth func- 
tion a{x, k) let the unique decomposition of a given by 

a(x, k) — a (x, k') + ai(x, k')k 3 + a 2 (x, k)(v \ k \ —u). 

Then we set Tj(a) = a, , i = 0, 1, 2. 

Similar results hold true for the other scalar semi- classical measures. 

Theorem 1.2. Calderon boundary condition case Using the same framework 
as in Theorem 1.1, but with Calderon boundary condition, the associated semi clas- 
sical measure jl ext (x,k), corresponding to the exterior part, is supported on the set 

U = {(x,k) , oj ext =uj} U {(x,k) , u ext = uj) . (1.19) 

Furthermore, it has the form 
' ji ext {x, k) = ^ e f-\x, k)b e f'\x, k) ® b extM {x, k) + ^ e f- 2 (x, k)b ext > 2 {x, k) 

^ 2 *{x, k) + li- t,l {x, k)b ext '\x, k) ® b extM (x, k) + n e * ta {x, k)b ext ' 2 (x, k) (1.20) 

®b e _ xt ' 2 *(x,k) 

t ext ,1 ext,2 , 7 • , ■ , j ,j , 

where u + ' , u + are two scalar positive measures supported on the set 
{(x,k) , u ext — lo], and y, ext ' , y, ext " , are two scalar positive measures supported 
on the set {(x,k) , uj ext = uj} . b^' 1 , fe^*' 2 (resp. tf**' 1 ,b ext ' 2 ) are two eigenvec- 
tors of the exterior dispersion matrix 

3 

L ext (x,k) = Y,( Aextfi y lk 3 A3 > 

3 = 1 

corresponding to the eigenvalue Lo ext (resp. uj ext ). 

The scalar transport equation, for the first scalar positive measure /i^*' 1 is given by 
Vfe^.V^T^-V^.V^r*' 1 = « ext fc3[<^ 1 ^3 =fer ,-+^ 1 5 fe3=fert , + ]4 3= o .(1.21) 

Above v e +*^ , v+*p~ are scalar positive measures associated with the semiclassical 

measure v ext corresponding to the boundary term u ext ' e ' ^(x 1 , 0). The wave vector 
k ext ^{k') = (k', kl x * ,=b ) is defined by 



3 v ' V v ext (x',0) 2 
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and v ext (x) = 



1 



is the propagation speed for the exterior problem. 



e ext (x)r] ext (x) 

Similar results hold true for the other scalar positive measures. 

Similarly, for the interior problem ( x 3 < 0), with the following interior dispersion 
matrix 



L mt (x,k) = J2(A int '°)- 1 k j A j 
i=i 

the corresponding semi classical measure jl mt {x,k), is supported on the set 

U={(x,k) , uj'f =lu} U {(x,k) , u™ 1 =uj} . (1.22) 
Furthermore, it has the form 
■ ji int {x,k) = fi™ t ' 1 (x,k)b™ t > 1 (x,k)®b™ t ' 1 *(x,k) + ^(x^^^ix^) 

(x, k) + M !_" M (x, k)b M '\x, k) ® &!_" M *(x, k) + ^(x, k)b M - 2 {x, k) (1.23) 
®b l ^ 2 *{x,k) 



i int. 2* / 



where ^i" 1 *' 1 , /i™*' 2 are two scalar positive measures supported on the set 
{(x,k) , uj" lt — ui} and fi™*' 1 ,fi mt ' 2 , are two scalar positive measures supported 



on the set |(x,fc) , u>_ = ui>. o + ,o + (resp. o_ ,b_ J are two eigen- 
vectors of the matrix L rnt (x,k) given above, corresponding to the eigenvalue lo™ 1 
(resp.uj™ 1 ). 

The scalar transport equation for the first positive measure fi 1 ™*' 1 is given by 



Y7 , Ant Y7 , int,l Y7 int ^7 ,mt 



Jul v" 1 v7 ,.,**»* V7, „*»»*> 1 V eXt ks\lf nt '' i ' 8 

intA intA 



1 int,l i 



7 mc,± tnz.i 7 - 1 j , j-i -7-7 

where v + ^ , v + A are scalar measures associated to the semiclassical measure 



"int 



„int,±{U\ 



corresponding to the boundary term u 
(k') — (k 1 , fcg"*^) is defined by 



int, £. 



<* fc ,= fc S»*.+ ]k,=o(1.24) 

•lassical measure 
(x',0), and the tangential vector 



1 int,± ( 1 



(x',0) = ± 



k' 2 



where v mt (x) — 



y v int (x',0) 2 
is the propagation speed for the interior problem. 



^fe mt (x)v mt {x) 

Similar results hold true for the other scalar positive measures. Finally, setting 

(1 o\ 

10 

10 
10 
\0 0/ 
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we have the following relation 

i) mt = Mi) ext . 

By adapting the proofs of the above two main theorems, we are also able to deal 
with the curved interface case. We sketch the proof at the end of Section III. 

2. Prerequesites on semi classical measures 

In this section, we recall some properties of semi classical measures which are useful 
in the analysis of high frequency propagation problems. For more details, we refer 
to [11], [15], [18], [20], [19]. 

Let / : R d — ► R n be in L 2 (R d ) n . Its (unsealed) semi classical transform is then 
defined as 

W(x,k) = -±- } [ e*y- k f(x-y/2)®r(x + y/2)dy. (2.1) 

\ Z7T ) JR d 

Its scalar semi classical transform is w(x, k) = Tr(W(x, k)). The function / can be 
scalar (n = 1), or vector- valued (/* denotes the transposed conjugated of the vector 
/). In the latter case its semi classical transform is an hermitian n x n matrix. 
We want to consider the semi classical transform of high frequency waves, i.e of 
functions f £ (x) which are oscillating on a given scale e, such that e — > 0. Our expo- 
sition follows the ideas of P. Gerard [25], [19]. Therefore, we consider the rescaled 
semi classical transform, at the scale e 

w£ ^ k ) = TTY d l e^ k f{x-ey/2)®f*{x + ey/2)dy. (2.2) 

\ Z7r ) JR d 



Proposition 2.1. Let the family f e be uniformly bounded in L 2 (R d ) n . Then, upon 
extracting a subsequence, the semi classical transform W £ converges weakly to a dis- 
tribution W(x, k) <G S (R d x R d ) n , such that Tr W(x, k) is a non-negative measure 
of bounded total mass (in the case n = d). 

Let a(x, k) be a test function in S(R d x R d ), where x € R d is the spatial variable, 
and k e R d is the momentum, or also the dual variable to x in Fourier space. Then 

<a,W £ >={a w {x,eD)f,f £ ) (2.3) 

where <,> is the duality product between S' (R d ) and S(R d ), (,) is the L 2 (R d ) 
inner product, and the Weyl operator a w (x, sD) is defined by 

[a^{x,eD)]f{x) = -^— d f , ek) f \y)e^ dkdy 

i r _u (2 ' 4) 

1 f ~ f x + y y-x 
(2tt£) J Rd 2 e 
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Here a is the Fourier transform of a(x, k) in the variable k only, 

a(x,y) = [ e- tk - y a(x,k)dk (2.5) 



and this operator is bounded on L 2 (R d ), uniformly in e, 

||a"'(x,eD)|| L 2 (R d ) _ +L 2 (R d ) < c(a) . (2.6) 
We also introduce the pseudo differential operator at the scale e, a(x,eD) by 

[a(x,eD)f}(x) = -±- } [ e tx - k a(x,ek)f(k)dk. (2.7) 

Again, one can show that the operators a(x, eD) are uniformly bounded on L 2 (R d ); 
there exists a constant c(a) > independent of e € (0, 1) (but depending on the 
function a) so that 

\\a(x,eD)\\ L 2 m ^ L 2 m < c(a) (2.8) 
and furthermore, it satisfies for any s > 

£ s ||a(a^eZ))|| i j-a (IR ci ) ^ I ,2 (K<1) < c s (a) (2.9) 

and 

e s \\a(x,eD)\\ L 2^)^ H s(Rd) < c s (a) . (2.10) 

The important point is that 

\\a(x,sD) - a w (x,eD)\\ L2{Rd) ^ L 2 {Rd) (2.11) 

as e — > 0, so that the two quantizations are asymptotically equivalent. 

With the above notations, one has the following link between pseudo differential 

theory and semi classical transforms 

lim (a{x,eD)f,f £ ) =< a,W >=Tr I a(x,k)W(dx,dk) (2.12) 

£ ► J 

(where we have also included the vectorial case). 

We shall also need the following results, from pseudo differential calculus (adapted 
at the scale e) 

Lemma 1. The product of two operators a(x,eD), b(x,eD) can be written as 

b(x, eD)a(x, eD) = (ba)(x, eD) + s/i(V k b.V x a)(x, eD) + e 2 Q e (2.13) 
where the operators Qe are uniformly bounded on L 2 with respect to e. 

Lemma 2. (Localisation) Let f £ (x) be a uniformly boundedfamiliy of functions 
in L 2 , and let ^if(x,k) be any limit semi classical measure. Let </)(x) be a smooth 
function. Then the semi classical measure of the family g £ (x) — <f>{x)f £ (x) is 
\(j)(x)\ 2 fj,f(x, k). Moreover, let f £ , g £ be two uniformly bounded families of L 2 func- 
tions which coincide in an open neighbourhood of a point x . Then any limit semi 
classical measure n / and /i g coincide in this neighbourhood. 
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3. Proofs of the Theorems 

3.1. Proof of Theorem (1.1), Perfect boundary condition case 

We consider the time harmonic form of Maxwell's system (1.11) or equivalently 
(1.6), in the half space of K 3 , (x 3 > 0), where E £ = (Ef,E e ,E§), and x = (x\x 3 ), 
x' G 

the outward unit normal vector, i.e. n = - 
equivalent to 

Ef = 



2 , x 3 > 0, with a perfect conductor boundary condition ft A E £ = 0, n being 



■k, which in our flat boundary case is 



Ef = 



(3.1) 



We set E £ to be zero in the lower half space x 3 < 0, and thus Maxwell's system 
(1.11) or (1.6) can be rewritten as 

3 



IW 

e 



A°(x)u £ + Y,A J 

3 = 1 



; dU £ 

dxj 



+ Cu £ = f{x) + A b u £ {x', 0) ® S X3=0 



(3.2) 



where the "boundary" matrix Ab is given by 

/0000 0\ 
0000 
0000 
0-10 
1 

\o o o o o oy 



A b 

In fact, let us recall that, 

rotE £ = V A E £ = 



(3.3) 




■d 2 E§-d 3 Ef 
d 3 E £ - d x E £ 
^El-dsEf, 



d 2 E £ 




'00 

-1 | d!E £ + 
1 



d,E £ . 



(3.4) 



As n A E £ — 0, we have that 



'd 2 E £ -d 3 E £ ' 
ml E V /. | 8 3 Ef - d!E £ 

d x E £ - d 3 E £ 



For the magnetic field H £ , let f2 C 
one has 



be a open domain dc . 



iotH £ .(pdx= / H £ .rotipdx— / (n A H £ ).tp dx 



(3.5) 

ThcnVy e C c °°(0) 3 , 
(3.6) 
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8> 5 X3=0 \ 


I- 




0^X3=0 






o / 



Extending H £ by zero in the full space, we have 

/ TotH £ .tpdx= H £ .vottpdx- (nAH £ ).tpdx. (3.7) 

JR3 J R 3 J K 3 

Here n = —k, thus we get 

. / d 2 H £ - d 3 H £ 

rotH £ = d 3 H £ - ftflf 

(3.8) 

fH £ ®5 X3 = \ (H £ \ 
rotH £ - I JTf ® ^ 3= o I - ™tii £ + I fff I ® <f Xs=0 

and using all the above notations, we get (3.2). 

Let 6{x) be a test function with compact support that is equal to one on a compact 
set K. We multiply u £ by 9(x), and thus we can define the semi classical measure 
jl on K for the family Qu £ , that we assume uniformly bounded in L 2 . 

More precisely, set 

u e ' 6 {x) = 6{x)u £ {x) (3.9) 

and let u £ > 9 (x',0) its boundary value, which is meaningfull in some negative 
Sobolev space, see [7] for instance. We shall assume that u £ ' 9 are uniformly bounded 
in L 2 (R 2 ) and that u e > 9 (x , 0)^3=0 are uniformly bounded in i?" 1 / 2 ""^) ( see 
[8] or [7]). 

We let (after having possibly extracted a suitable sub-sequence) jl and v be the 
(matrix valued) semi classical measures of u £ ' " and u e ' (x 1 , 0) rcsp. 
Now Maxwell system can be rewritten, with the cuttof function 6, as 

" -,du £ ' 6 .< dO 



iluA°(x)u £ > 9 + e V A> ^— - e V A i -^-u e (x) + eC{x)u £ ' 9 

OXj OXj 



(3.10) 



Let a(x, fc) be a matrix-valued test function with compact support in K, with 
respect to x. Applying the operator as = a(x,sD) on both sides of (3.10), and 
taking the inner product with u £ > 9 , wc get 



(a £ [iu>A°(x)u e > 9 +eYA3 -^zl - £ y^ ^-u £ (x) + eC(x)u £ > 9 ] , u £ ' 9 ) 
z — ' oxj z — ' 



9xj pi d Xj (3.H) 



= e(as{f> °{x%u £ ' 9 ) + e(a e [A b u £ ' e (x',0) ® 5 Xa=0 ],u e > 9 ) . 

This is well defined in view of the usual rules of pscudo differential calculus (at the 
scale e). 
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To evaluate the limit of the second term of the right hand side in (3.11), let us set 

v £ := A b u e > 6 (x ,O)0^3 = o 

and thus we get that v £ is uniformly bounded in _ff _1 / 2_0! for any a > 0, with 
a < 1/2. 
Next, note that 



\e(a £ {v%u £ > e )\ < e\a £ [v £ }\ | |u e ' °\ \ L * m < 

e\\a £ \\H->— -*L2\\v £ \\H-A\ u£,d \\L2(W3) < C^|||w £ ||iJ-3||M e '^||i,2 (K 3) . 



(3.12) 



Thus if we choose s = —1/2 — a, we get 



1/2 



\\£(a e [v £ },u £ > d )\\ L 2 {R3) <^ , a>0. 
It follows that choosing a < |, one has 

lim e(a £ [v £ ],u £ > 6 ) = . 

£ ► 



(3.13) 



Let us also show that the other terms in (3.11) are bounded (uniformly with respect 
to e). Indeed, for the third term of (3.11), one has 



\(ae[-e^Ai^u%u e > 9 )\ < e\a £ [J2 ^ -^u £ }\ \\u £ ^\\ L2{R3) < 

3 = 1 2 j=l 3 



(3.14) 



09 



e||ae||L2 (R3) || ^ A 3 ^-M e || L 2 (R3) ||M e ' e || L 2 (R3) < ec 



3 = 1 



dx 



and thus 



£ lim (-ea £ [E A> -j£-u%u s > °) = 



i=i 



(3.15) 



and similary for the terms e{a £ [C{x)u £ ' ], u £ ' ), and e(a e [/ £ ' ], u e ' ). 
For the second term on the left hand side of (3.11), we set 

a 2 :=e(a £ [±A^},u £ ^, 

3 = 1 3 

a £ = a(x, eD) , 



(3.16) 
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and using the product rule (2.13), we get that 

' a 2 = (a £ [b £ u £ > 9 ,u £ > 9 ) = (a(x,eD)(b(x,eD)[u £ > d ],u £ > e ) w {{a £ b e )[u £ '% u £ ' 9 ) 

3 

Tr / / a(x,k)b(x,k)jl{dx,dk) = -Tr / / a(x, k){i kjA j )jl(dx, dk). 

JR3 J Ml JR3 J Ml ~i 

x k x k J— 1 

For the first term of (3.11), we set 

' a x := {a e [iujA (x)u £ ' e ],u £ > e ) , 

a £ = a(x,sD), (3.18) 
b £ = A a {x) 
and thus, one has 

oi := (ae [iwb(x, eD)u £ ' 9 },u £ > 9 )^ iw(a £ b £ [u e ' 9 ] , u £ ' ) 



(3.17) 



iwTr 



o(x, k)b(x, k)fl(dx, dk) = iujTt 



a(x, k)A°(x)fl(dx, dk). 



(3.19) 



Thus all in all, passing to the limit in (3.11), we get 



Tr / I a[ 

[3 ./»3 



(3.20) 



i(x, k)[iujA° - i V" kjA j ]fl(dx, dk) = 

for all matrix valued test function a(x, fc) £ <S(R 3 x R 3 ), which is equivalent to (in 
S>) 

3 

{iwA° -iJ2kjA j )fi = (3.21) 



or 



Let us set 



(-A°)[i (A°)- 1 k j A j -iujld]fl = 
3 = 1 



L(x,k) = ( A °r lk J AJ ■ 

3 = 1 



(3.22) 



(3.23) 



In order to find the eigenvectors of the 6x6 matrix L(x,k), we shall use an or- 
thonormal propagation basis of R 3 . We denote by (k, z 1 (k), z 2 (kj) the orthonormal 
propagation triple consisting of the direction of propagation k = kj | k | and two 
transverse unit vectors z 1 (k) 1 z 2 {k). In polar coordinates, they are, see for more 
details [19], [20] 

' sin 6 cos d> \ / cos 9 cos <t>\ I — sin < 

z 2 (k) = ( cos0 | ru-n 



k = -j-^-j- = I sin sin J , z 1 (fc) = I cos (7 sin ( 
cos# / V — sin( 



I 



I 
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where |fc| = {k 2 + k\ + fcf) 1 / 2 . 

Then the eigenvectors (which belong to M 6 ) of the matrix L(x, k) are given by [see 
[19], [20]] 



6j = ^=(fe,0), 6§ = -L(0,fc) 



l V2^ 'V^ 



b 2 - (—z 2 



s 1 ) 



(3.25) 



bl = 



(-/=** > 7= 



z 2) b 2 ={ '2' 



The eigenvectors 6j and 6q represent the non-propagating longitudinal and the other 
eigenvectors correspond to transverse modes of propagation with respect the speed 
of propagation v. These eigenvectors correspond to the eigenvalues listed in the 
following Lemma, whose proof follows from [20], [19]. 



Lemma 3. The semi classical measure jl is supported on the set (recall that we 
assume that the frequency uj ^ 0) 

U = {(x, k) , uj + = uj} U {(x, k) , w_ = uj} (3.26) 

where v(x) = — = is the propagation speed, and ujq = ujo(x, k) = , uj + = 

uj+(x,k) = v(x)\k\ , uj- = oj-(x,k) = — v(x)\k\ are the eigenvalues (of constant 
multiplicity two) of the dispersion matrix L. 

II follows that the semi classical measure fl(x, k) has the form 

jl(x, k) = n\(x, k)b\(x, k) <g> b^(x, k) + n%(x, k)b 2 (x, k) <g> b 2 ?(x, k) 

(3.27) 

+fj}_(x, k)b x _{x, k) <g> feL*(x, k) + fi 2 _(x, k)b 2 _(x, k) <g> bt*(x, k) 

where n\ are two scalar positive measures supported on the set 

{(x, k) , uj + = uj}, and fi\_ are two scalar positive measures supported on the 
set {(x, k) , uj- = uj}. b\ ,b\ (resp. b\_ b 2 _) are the two eigenvectors of the matrix 
L(x, k) given by (3.25), corresponding to the eigenvalue uj+ (resp. uj-). 
In view of the above reduction, we are led to find the transport equations for each 
of these four scalar semi classical measures. 

For this purpose, using the equation (3.10), we have the following identity 
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= ioj{a e [u £ ' 6 ] , it 6 ' 9 ) - iiv{a e kt £ ' 6 ] , u £ ' °) = (iuja £ [u e ' 6 ] , u £ ' 6 ) + (a £ [u £ ' °] , iwu £ ' 9 ) 
= (a £ [-e £ (^-i^i^! + e £ (A^{x)A^u £ {x) 



3=1 J 3=1 

-e(AY\x)C(x)u £ ^ 9 + e(A°)-\x)f^(x) + e(A°y 1 (x)A b u £ ^(x' ,0) (g> 6 X3=a ],u £ ^ 9 ) 



+ {ae[u £ A-s (A a )- 1 (x)A3?^—+e £ {A °)-\x)A^u £ (x) 
-e^ )- 1 ^)^^)^' 9 + e(A°)- 1 (x)f £ ^(x) + e{A°)-\x)A b u £ ' 9 (x\ 0) ® 4 3 =o) 



(3.28) 



= - £ (a e [ ^ (A°)~ 1 (x)A j ^—},u £ > 9 ) +e( £ — [(^W^' 



£,6»i u e,6K 



i = i 



i = i 



+ £ (a e [ ^ {A a )- l {x)A^u £ {x)] 1 u £ ) - s(ae[(A )- 1 (x)A b u £ ^(x, 0) ® 4 3 =o],« £ ' ) 



.7 = 1 



9a; 



90 



+e(ae[u £ > % ]T (A )- 1 ^ '^ U £ ) - e(a e [u £ ^ % (A )' 1 (x)C(x)u £ > %u £ > 6 ) 



3 = 1 



dx 



+e(a £ [u £ > % (A°y 1 (x)A b u £ ^ 9 (x , 0) ® <5 X3=0 ) . 



Recalling that the function 9 is equals to one identically on the support of a(x, k), 
the third, fifth and sixth terms vanish at the limit, and thus the last equation can 
be rewritten as 

3 o fO 3 



e(a £ [ £ { A°)-\x)A^^—lu^ e )+e{ £ Ap")- 1 ^^^^) 

x 3 . = i 3 



3 = 1 



= s(a £ [(A )- 1 (x)A b u £ ^(x, 0) ® 5 Xs=0 ],« e ' ) 
+ £ (a e [ u £ ' (A°)- 1 (.T)A bU £ ^( a; ', 0) ® 5 XS=0 ) . 



(3.29) 



Using the product rule (2.13), it follows that 

3=1 



dxj dxj 
j J=1 j 



where 4>q ,<j>\, are given by 

3 3 

cj> a (x,k) =ia{x,k)^2{A°)- 1 {x)k J A J - i^kjA 3 {A°)- 1 {x)a{x,k) , 



3 = 1 



3=1 

3 ~, -n\-i 3 



9a d(A°)-i Aju ^ Ai 0{A°)-\ ^ AUA ^ da 



(3.31) 



j,m=l 



dk m dx r ' 



-A j kj -J2 A3 ' 



3 = 1 



dxi 



3 = 1 
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and the operators Rs are uniformly bounded on L 2 . 

On one hand, using (3.30), and the two relations (2.9), (2.10), we pass to the limit 
in (3.29), as e — > 0, and obtain 

Tr / / <j> (x,k)a{x,k)ii{dx,dk) = 0, \/a(x,k) (3.32) 
which is already a known result (localization principle). 

On the other hand, dividing (3.29) by e, and passing to the limit as e — > 0, we get 

Tr / / <j) 1 (x,k)il{dx,dk)+ lim -(<j> {x, eD)u 6 ' °, u e > °) = lim B £ (a) (3.33) 
JR3 J Ml e — > £ £ — > 

where 

B e (o) = (a e [(A°)- 1 (x)i4 4 u e ' e ®5 I , = o],« e ' e )+(a e [« e ' e ],(A )- 1 (x)i4 4 u e ' fl (x')®J x , =0 ) .(3.34) 

In order to find the transport equation for the semi classical scalar measure jl, we 
use the orthonormal propagation basis, and we consider first a test function a(x, k) 
of the form 

a(x, k) = a + (x, k)d\{x, k) <g> d^(x, k) (3.35) 

where a + (x, k) is any scalar smooth function, and 

d\(x, k) = A°(x)b 1 + (x, k). (3.36) 

Recalling that A is a symmetric matrix, with the choice (3.35), we note then that 
0o vanishes, while <p\ becomes 

M x, k) = ± ^g^d\ (x,k)® d\* (x, *) ± A> kj 

m— 1 j — 1 

£ da +^ k ) A1{A o r i d i +{Xi k) d i r(a . ) fc) 



9xj 

m=l m j'=l 

m— 1 — 1 

3 a Jl 3 rW 1 * 

^^(A )- 1 ^ ® 4*(x, fc) - ^ A^A )-H\{x,k) ® ^-±-} - U + 12 + 13 . 



(3.37) 



l 



l 
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where 



E 

m— 1 



da + (x,k) x 



dk, n 



d\{x,k) ®d+(x,k) 



E ' 



E 



da+(a;, fc) 



M.2 = -7 — — — '-A j (A°)~ 1 d 1 i _(x, k) <g> c^*(x, fe) , 



»13 



1 dd\*{x,k)d{A )- 1 ' 



>E 

m— 1 

3 



dk„ 



dx r ' 



(3.38) 



- E ^^"H ® fc ) - E A\A^d\{x, *) ® ^} ■ 

We shall use the eigenvectors in the orthonormal basis (3.24), and the following 
normalization relations, 



' {A°b a ,bp) = 8 aP , 

< 

^ {Aib+,b + ) = vkj. 
We can then evaluate the first term 0n in (3.37). Indeed, we have 



(3.39) 



< 0n, n >=< 



da + (x,k), ,^^d(A°)-i 



dk n 



d\[x,k)®d+{x,k)Y^ v J. A J 'fc J -,^(a ; ,fc)6j F (a;,fc)0(b + (a;,fc)) 1 *) > .(3.40) 



Since 



d 



3 



■(E(^ r 1 M J <)=E 



It follows that 

3 



a(A°)- 

9x m 



■M^ + D^r'M^- (3-4i) 



Eipr'M^-B^M^^^^i'^ ■ (3.42) 



i=i 



i=i 



Using the eigenvectors of the dispersion matrix in the orthonormal basis (3.24), one 
has 



dx m 



(3.43) 



l 



l 
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Thus the first term in (3.37), becomes 

3 



(3.44) 



Using (3.39), it follows that 

(|^, tf) = (^AX.&i) = ^ ,m = 1,2,3 

(^+^r^+) = <"+( ^^rM) = o ,m = 1,2,3 (345) 



All in all, the first term of (3.37), becomes 

3 



m— 1 



For the second term 12 in (3.37), one has 

< 0i2, A >= - < 9a+ Q^ k) (A°)- l A>d\{x, k) ® fc), > 

3 J=1 ^ (3-47) 

= -<E ^^(buim^r^dui),^ > . 

j = l 3 

Using the eigenvectors of the dispersion matrix in the orthonormal basis (3.24), and 
(3.39), one has 



{b\,d\) = {b\,A%\) = l 

duo + 



((A°)-'A^d\,b\) = ((A )^A^A°b\,b\) = {m\,b\) H 



(3.48) 



and the second term </>i 2 in (3.37), becomes 

<^,,>=-<J2^^,,l> • (3.49) 
j=i ■> 

For the third term ^13 in (3.37), we shall show that 

<0i3,/*>=O. (3.50) 



l 



l 
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We need to deal with the following term 

m— 1 m j — 1 

™ =1 3=1 (3.51) 

-E AJ '-Sd~ d + (iB ' fc) ® d + (iB ' fc) 

- ^^(A )- 1 —^- ® d + ( x ' fc ) - E A- 7 '(A°) _1 d]j.(ar, fc) ® -^}, ^(x, fc)4(s, fc) <8> d^fo fc) > . 



Set 

m=l " j'=l m=l j'=l m 

-D^^^^^^^-D^.^^ )- 1 ^)^^) (3-52) 
i=i j'=i 

We can rewrite (3.51) as 

< ^i 3 ,M>-<a+[T],M+ > • (3.53) 
For the first term in (3.52), we use (3.24) and (3.39), to get 

m—1 ' j—1 m—1 j — 1 

m—1 m j—1 

= E «.E^^O = 

TO=1 m J = l 

and thus the first term in (3.52) vanishes. For the last term in (3.52), we use again 
(3.39) to get 

' E(^.^(^4)(^.^) = E(6i,^(A r^X)(6i.^) 



i=i j=i 



E(^.^ + )(^. ^) = E^ & +< |J) = ^) (3-55) 

3=1 ' 3 3 = 1 3 3 

duj + , <9di s 



X-^±( b i 



i 



i 
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Thus (3.52) becomes 

j — 1 711 — 1 J = l J= 1 

For the last term in (3.56), we use (3.24) and (3.39) to get 

3 du+dd^^ _ ^ dw+ fhl ^ _ _^ dw+ fhl d [A°b\] ^ 

(3.57) 



a* j - 2- 9fc, j _ ^ dkj [b+ ' d& ' 

3 = 1 3 = 1 J=l 



-t^{(1.^) + «.^)}-t-^»l.^»- 

For the second and third terms in (3.56), we use (3.24) and (3.39) to get 

- W ,£{^4 + ^vg}>-<*.£*^> 

3=i ^ 1 ' j=i 1 (3.58) 



For the first term in (3.56), we use (3.24) and (3.39) to get 

(+ 'fe^i M «^ ^~ { +, h h d * m dk - ] , 3 59) 

j = l m=l 



But 



Thus (3.59) becomes 

3 3 „, otl 3 3 



v+'Z^Z^ ftr™ dfc„/ v+ 'f-f ^ dx m dk m 



J — 1 777—1 J — 1 777—1 



j — 1 ra—\ ' m — 1 



Thus all in all, (3.56) becomes 

— «.E^S»-(^.§. + ^ + ,,»§.. (, 62) 

3=1 



I 



I 
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For the second and third terms in (3.62), using the fact that b\ an eigenvector of 
the dispersion matrix, one has 

° 3 j=i 3 3 = 1 3=1 0ltj (3.63) 

d . <9w +ll db\ 

= dk- {uJ+b+) = ^k- b++UJ+ dk; 
which implies that 

9kj dkj dkj 

and thus second and third terms in (3.62) become 

-(^■§)^Ci.^>-(1.t^f-^Vf). (3.4, 
Thus all in all, we have 

3 = 1 J 

which yields (3.50). 

Now, using (3.46), (3.49), (3.50), and integrating by parts, (3.33) becomes 

3 3 



(3.65) 



< <Au,M > + < 0ia.li > + < 013,/* >= E < ft^^.Mf > -L < afc7^'M + > (3 



=< a + , VfeW+.V^^I - V x w + .Vfc/ii >= lim i? £ (a) . 

s — > 

There remains to determine the right hand side of (3.66). Recall first that 

B £ (a) = {a £ [(A' ) )- 1 {x)A b u^ d (x', 0)®J X3=0 ], u £ > [« £ ' % (A°)~ 1 (x)A b u £ > 6 (x, 0)®<S X3=0 ).(3.67) 

Note that each term in (3.67) is of order £ _1 / 2_Q for any a > as can be seen 
from the H s estimates, since u £ >®(x ,0) <g> £ X3 =o is uniformly bounded in H s for 
s = — 1/2 — a, for any a > 0. 

To get the limit of (3.67), we shall first use a special class of matrices a(x, k) of the 
form 

a(x, k) = a(x, k)[L(x, k) — u>I] (3.68) 
where L(x, k) is the dispersion matrix (3.23) and for any matrix a(x, k) satisfying 
a(x, k) [L(x, k) - luI] = [L* (x, k) - wl]a(x, k) . (3.69) 



I 



I 
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Using the test function (3.68) and the product rule (2.13), the first term of (3.67) 
can be worked as follows 

(a e [(A°)- 1 (x)A b u s ' 9 ® 5 X3=0 },u £ > 9 ) 
= ^([L*(x, k) — ujl]a(x, eD) [(A°)- 1 (x)A b u £ ' d ®S X3=0 ] ," £ ' ) 
~ (([L*(x,k) -ujI}(x,eD)&(x,eD)[(A°)- 1 (x)A b u £ ^ ® 8 X3=0 ], u e > ) 
- t ((V fc [L* (x,k)- u>I] . V x 5) (x, eD) [(A )- 1 (x) A b u £ > 9 ® S X3=0 ] , M e ' + e 2 Q 

- ({a{x, eD)[(A°)- 1 (x)A b u £ ' 9 ® 4 3 =o], [£(*> fc) - w/](a;, £ £>)u £ ' 9 ) 
-| ((V fc [L*(z, fc) - wil.VxfiJCa:, e-D)[( J 4°) _1 (a;)A(,u £ ' ® 5 X3=0 ],« e ' 9 ) + e 2 Q £ 



(3.70) 



with a term Q e uniformly bounded. 

The two last terms of the above formulae are uniformly bounded and vanishes to 
the limit. 

Indeed, for the first term, recall that 



^ dxi 



(3.71) 



j=l ~ 3 

We then use (3.11) and (3.71), to rewrite (3.70) as 

(a s [(A°)- 1 (x)A b u^ 9 ® <5 X3=0 ],« £ ' ) 
(~a(x, eD)((A a )- 1 (x)A b u £ ^ 9 (x) ® 5 X3=0 ), -(A°) _1 (x)/ e ' ^(a;) + -(A°) _1 (a;)A b u e ' V 



i=i 



(a(^,£ J D)((^ )- 1 (x)A 6 u e ' w (x ) ® 4 3=0 ), +-(^°)- 1 (x)A fcU £ ' p (x ) ® <5 X3=0 

2 



For the second term of (3.67), still using the test function (3.68) and the product 
rule (2.13), we get that 



(a £ [u £ i °], (A°)- 1 (x)A b u £ ^(x) ® S X3=a ) 
= ({h[L{x,k) - LoI]{x,eD)[u £ ^ 9 ], \ i (A°)- 1 (x)A b u £ ^ 9 
~ (ia{x,eD) [[L(x, k) - ul]u e £ ] , ^(A a )~ 1 (x)A b u £ > 9 ® <5 : 
-| ((V fe a.V x [-L(z, fc) - o;7])(a;, eD), ^(A )~ 1 (x)A b u £ > 9 ® <5 X3=0 ] ) 



2:3=0 
23=0 



) 



(3.73) 



e 2 R £ . 



2:3=0 



(3.72) 
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with Rs uniformly bounded. 

We use (3.11) and (3.71) to rewrite (3.73) as 



' (a e [u E > ],(A°)- 1 (x)Ai,u £ > {x)®5 Xa=o ) = (a(x, eD)( £ -(A )-\x)A b u £ > d (x ) ® S X3=0 



+ i(A r 1 (x)A h ^-u £ ^, [(A a )-\x)A b u £ > d (x)®S X3=0 



) 



Using these asymptotic expansions, passing to limit in (3.70), as e — > 0, we obtain 
finally 



lim B £ (a) = . 

£ — > 



(3.75) 



Now, we consider the general case of test functions in order to pass to the limit in 
the boundary term. For this purpose, we note that it is possible to write every test 
function a+ as 



a + (x, k) = a (x, k ) + a\(x,k )k 3 + 122(2;, fc)(u|fc| — lo) 



(3.76) 



where k = (k ,k 3 ) and ao , a\ and 03 are scalar test functions, uniquely determined 
by a+. For this point, we refer to [23]. 
In view of (3.76), we shall set 



T (a+) = ao, Ti(a+) = ai and T 2 (a+) = a 2 . 
Then any a of the form (3.35) can be written as 

a(x, k) = (a (x, k ) + a\{x,k )k3)A°(x) 



(3.77) 



-(a 2 (x, k)d\ ® d 1 ^ 



a (x,k ) + ai(x,k )fc 3 ^ 1 
v\k\ + uo 



(3.78) 



+ a (x,k) + a 1 (x,k)k 3 _ 4 ^ _ ^ k) 

Now, we note that the spectral representation of the matrix L — uoI can be written 

as 



L 



jI = (gj+ - oj)b\ ® rf^* + (oj- - uj'jb^ ® rfL* + (w - w)6° ® d° 



(3.79) 



Recall that that the last term in (3.78) has the same form of the test function 
a = d[L — u)I] of (3.68), and thus we can conclude for the limit of this term and we 
find that 



lim B e {a 2 ) = . 

£ — > 



(3.80) 



Therefore, it is enough to find the limits for B s (a) only for the first two terms. 
For this purpose, denote by a = a (x, k )A°(x) the first term in (3.78). Multiplying 
it by a suitable cutoff function, (f>(e 3 k 3 ), with support compact, equal to one on a 
neighbourhood of zero, set 

a'e = a'(f>{e 3 k 3 ) = [a {x, k')A°(x)}<j)(£ 3 k 3 ) . (3.81) 



(3.74) 



l 



l 
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Using the product rule, the first term B\e{a ) leads to 

B 1E (a") = B £ (a £ [{A°)^A b u £ ® 8 Xs=o ],u E > ) 

~ (a (x,sD)A (x)(l>(e 3 k 3 )[(A )- 1 A b u^ d ®6 X3=0 ],u £ ' 9 ) 



(3.82) 



u e ' e *{x)dx 



dk 



e ik - x a {x, ek')<j)(e 3 k 3 )A b u^ (k') 



(2tt) 3 

and similarly for the second term 

B 2 e(a") =B e {a £ [u £ ^ ], (A°)~ 1 A b u £ > {x) <g> £> X3 =o]) 



u £ >"*(x )dx 



dk 



■e ik> - x ' a {x ,0, /c')0(e 3 fe)Aw e ' 6 {k) 



(3.83) 



(2tt) 3 

Thus for the first term in (3.67), and for the test function written as (3.76), one 
has 



lim B £ {a)= lim B 1£ (a ) + lim B 2 e{a ) = Tr / A b a (x , 0, k )dv (3.84) 
e^O e — > e^O J 

where v is the semi classical measure of the boundary term u £ (x , 0) . 

For the second term in (3.67), denoting by a = ao(x, k )ksA°(x), in the same way, 

we have 

B 1£ (a") = B £ (a'^A )- 1 A b u £ ' ® 5 },u £ > ) 
= (ai(x, e^)A o ( 2 ;)0( e 3 fc 3 )[(A o )- 1 A 6U e ' ® 4 3 = ], « e ' ) 



(3.85) 



Also 



B 2£ (a")=B £ (%[w e ^], {A^A h u £ ^{x) ® 5 I8=0 ]) 
dfo £ £)fc 3 A> £ > (yl )- 1 ^' V) ® 5 X3=0 



(3.86) 



= y u £ > d *(x )dx J -^e ik ' ai (x ,O,k')<t>(e 3 k 3 )A bU z,0(k) . 
Thus, we have 



Bie(a ) +B 2e (a' ) 



(2tt)3 
: / du £ > * 



+ lj u £ > 9 *(x)dx 
Passing to limit in (3.87), we get 



dx 3 

dk 
(2^f 



/ %x^ {x)dX J ^ elk - Xa i( x > £k 'W £3k 3) A ^( k ') 



> i f)u £ i 

»J k x ai (x ,0,k )0(£ 3 fc 3 )A^^(fc) . 



(3.87) 



r 2 

lim (B 1£ (a") + B 2 e(a"j) = -Tr / [V^ kjA 3 — ljA°(x ,0)]ai(x ,0, k')di> . (3.88) 
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Thus all in all, we get the limit of the boundary term (3.67) as 

f 2 
lim B £ (a) = Tr / [A b (x , 0, k')a (x', 0, k')-(Y] kjA j -ljA°(x , 0))ai(a;', 0, k')]di) .(3.89) 
e - J 

Note that if the test function a(x, k) inside is supported away from x 3 = this 
limit equals zero, as it should be. 

Now, because we are using special test functions satisfying (3.35), and since we are 
dealing only with we can as well assume that we are only seing the following 
part of v given by 



v ~ vl+b\{k + ) ® b 1 + *(k+) + ul fj+ b\{k+) ® b\*(k-) + 

+vi a+ b\{k-) ® bi*(k + ) + 4 + b\{k-) ® bi*(k-). 

This follows from the corresponding localization principle on the boundary. 
Next, note that 



J2(A°)- 1 k j A j - ujld 

3 = 1 



for scalar measures. 
We have also 



(A b b 1 + (k),b\(k)) = 0, 
(A 3 b\(k),b 1 + (k)) = vkz , 
(A 3 b\(k+),b\(k-)) = 0, 
(A b b\(k+),b 1 + (k-)) = , 

2 

kjA j - wA )^(fc ± ) = -fc±A 3 ^(fc ± ) 



(3.90) 



Using (3.90), (3.78) and (3.76), the term (3.89) becomes 

| lim^_Be(a)= J is^+(dx ,dk )vk^ a\ + J Vp + (dx , dk )vk^k^a\ (3.91) 



Recall that v(x) 



k' e K 2 , and the wave vector fc ± (fc') = (k' , kf) is defined by 



y/e(x)r)(x) 
tor k'- 

kf(x',0) = ± 



is the propagation speed, the tangential vector 



v(x\oy 



-k' 2 



i 



i 
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By using formulas linked with the wave vectors, and in particuliar definitions given 
in (3.77), the above formulae reduces to 

lim B £ (a) = 
e — > , 

, , r , , ( 3 - 92 ) 

vi + (dx ,dk )vk 3 k 3 T 1 (a+)(x',0,k') + / v} + (dx ,dk )vkt,k+Ti{a + ){x' , 0, k'). 

Combining (3.92) and (3.66), we get the following distributional form of the trans- 
port equation for the (scalar) positive measure v+(x, k) 

VfcW+.V x/ u+ - V x uj + .V k ^]_ = wfc 3 fc 3 [^ + Ti(5 fe3=fc - + v}+Ti8 k3=k +]d X3 = . (3.93) 

The other semi-classical measures in the formula are also dealt with in the same 
way as above, and we get 

VfcW+.V x /U+ - V x LJ+SJ k ^ 2 + = vk 3 k 3 [vl + Ti5 k3=k - + Vp + Ti6 k3=k +]6 X3 =o (3.94) 
Vtw+.V^L - V x uj + .VkH- = vhk^v^TxS^,,- + Vp-Tid k3=k +]8 X3 = (3.95) 
VfeW+.VxM^ - VzUJ+.VkH- = vkak3[Va-Ti5 k3=k - + t>p_Tid k3=k +]5 X3 = (3.96) 

3.2. Proof of Theorem 1.2, Calderon type boundary condition 

In this case, for the exterior problem (1.13) (given in R+), extending by zero in the 
full space R 3 , we have the following eikonal equation for the exterior problem 



iwA ext '°(x){u ext ' £ ' e ) + sJ2A j {U ' } - e J^A i jf-u ext ' E {x) 

oxj oxj 



= 1 J J=l 



(3.97) 



+eC ext (x)(u ext < e ' 6 ) = ef> 9 (x) + eA 3 u ext - £ ' (x',O) <g> 5 X3 = - 



Note that, on the contrary of the perfect conductor case, we have not at this level 
taken into account Calderon transmission condition. We have also includede in the 
exterior field the incident one, using the same notations. Above the matrix A ext '°(x) 
is given by 

^'"=( £ 7V°id) <398) 

where e ext ,rj ext , are smoth functions in C 1 (M 3 ), and the matrices A^ are given by 
(1.8), and the matrix C ext is given by 

T ext 



with <r ext a smooth function in C 1 (M 3 ), and u mt < 6 >® (x ,0) is the boundary term 
for the interior problem (i.e. x 3 < 0). In this case, we obtain that the dispersion 
matrix for the exterior problem is given by 

3 

L ext (x, k) = Y,(^ ext '°)~ 1 kjA j . (3.100) 



C EXt = ( ^ ° ) (3-99) 



l 



l 
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Recall that the matrix L ext has also three eigenvalues which constant multiplicity 
two. They are 



ext 




,w«f* = v ext \k'\ ,wff* = -v ext \k'\ 



where v ext (x) — -^^^=^^^= is the propagation speed for the exterior problem. 

\J e ext {x)rf- xt (x) 

As in the perfect conductor case, it follows that the associated semi classical measure 
fl ext (x, k) has the form 

{'li ext {x, k) — p? xt ' (x, k)b ext ' (x, k) ® b ext ' *{x, fc) + fi^*' (x, k)b ext '^{x, k) ® &^*' 2 *(x, fc) 
+ M l xt ' 1 (x,fc)&! xt ' 1 (x,A0®&^ 

where /z^f 4,1 , ^^*' 2 are two scalar positive measures supported on the set 
{(x, k) , w^?' = w}, and /il^*' 1 , /il xt ' 2 , are two scalar positive measures supported 
on the set {(x, k) , wf_ x * = w}. b ex ' ,b ext ' (resp. 61^*' ,b ext ' ) are the two eigen- 
vectors of the matrix L ext (x,k) given by (3.100), corresponding to the eigenvalue 
io ext (resp. lu^). 

The semi classical measure jl ext is supported on the set 

U = {(x,k) , oj ext = oj}u {(x,k) , u ext = oj) . (3.102) 
For instance, the transport equation for the first scalar measure is given by 
Vfc^V^^-V^'.VfcMf - 1 = v ex %Kf % 3=kr ,-+^ 1 S k3=kr . + ]S Xa=0 (3.103) 

where v ex + x ; t/ l g+' 1 are scalar measures corresponding to the boundary term 
u ext ' £ > 9 {x',0), and the wave vector k ext ^{k') = (k' ,k^ xt ' ± ) is defined by 



k ext ^(x',Q) = ±J—£—— - fc'2 . 

In fact, as in the previous sub-section, these (measures) coefficients come the fol- 
lowing decomposition of the boundary semiclassical measure (seing only the first 
part of the set U) 

~ext _ I/ «rt,l 6 e*M( jfe «rt,+ ) h ^*{k ext - + ) + U e a X ^b e f'\k eXt ' + ) ® b"f A *(k ext '-) + 

+v e £+b e f- 1 (k ext '-) ® &; xt ' 1 *(fc ex *'+) + u^'Hf' 1 ^"*'-) ® b ext ' U {k ext '-). 

For the interior problem (1.14), (given in K 3 .), we have the following eikonal equa- 
tion 



Q( u int,S,6} 3 gg 

£J (3.104) 



z W (A m *' )(x) («<"*■ £ ' e ) + £ V A J ( " ' ' ' -ej"A* ^-u mt - £ (x) 



+eC mt (x){u mt ' £ >°) =ef>°{x)+ eAl nt u ext - £ >° {x , 0)®5 X3=0 ■ 



l 



l 
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where the matrix A '(x) is given by 

where e mt 7 rf nt are smooth functions in C 1 (R 3 ), and the matrices are given by 
(1.8) and the matrix C mt is given by 

with a mt a smooth function in C 1 (M 3 ). 

In this case, the interior dispersion matrix is given by 

3 

L mt (x, k) = Y,(A int '°)- 1 k j A j (3.107) 

has three eigenvalues with constant multiplicity two. They are 
4 nt = ,Lu+=v mt \k'\ ,w_ = -v mt \k'\ 

where v mt (x) = — is the propagation speed of the interior problem. 

y/e mt {x)rf nt {x) 

Again, it follows that the associated semi classical measure ji mt (x, k) has the form 

f fl mt (x, k) = tff'Ux, Qbf'Hx, k) ® b™ tM {x, k) + nf^ix, fc)6" l *' 2 (x, k) ® &^^ 2 *(a;, jfe) 

<^ (3.108) 

1 V- M (*> fc)&!_" M (z, k) ® &!_" M >, fc) + M-*' 2 (z, k)b M - 2 {x, k) ® &!_" M *(x, jfe) 

where /i™*' 1 , /i™*' 2 are two scalar positive measures supported on the set 
{(x, fc) , = w}, and /x™*' 1 , /x!! 1 *' 2 , are two scalar positive measures supported 
on the set {(x, /c) , w_ = u>j. b + , o + (resp. b_ ,o_ ) are the two eigen- 
vectors of the matrix L mt (x,k) given by (3.107), corresponding to the eigenvalue 
w™' (resp. o;!! t *)s. 

The semi classical measure jl mt is supported on the set 

U = {(x,k) , w; n * =w} U , w mt = . (3.109) 

As an example, the transport equation for the first scalar measure is then 

Vfc^.Vxi/^-V^.Vfci/^* - V <nt fe[i/^ 1 5 ta=fc j»*.-+i/^' 1 5 fc3=fc i- t . + ]5 a! ,=o -(3.110) 

where v™^ 1 ,v™^ x are scalar measures associated with the semiclassical mea- 
sure corresponding to the boundary term u mt ' e ' ^(x 1 , 0), and the wave vector 
k mt ^{k') = (fc',fc™*^) is defined by 



Let us end by making some remarks about the scalar measures appearing the right 
hand sides of each transport equation, in the exterior as well as in the interior case. 



I 



I 
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Due to the Caldcron transmission condition, it follows that one has the following, 
on the boundary x = 

u int, £,6 — m u ext -. £ : O 



where 



M = 



It follows that (M 2 = M) 



(\ o o\ 
10 

10 
10 

\0 0/ 



i) int = Mv ext . 



For instance, to get the scalar measure v™+ x , it is enough to take the trace of the 
above relation with ci™*' 1 (/c m *' + ) ® d™*' 1 *(/c m * ,+ ) (where we use left eigenvectors) 
and we get in this way 



Tr(MD ext .d i ^ t ' 1 (k int '+) ® d l ^' 1 *(k mt '+)). 



■int,l* /> int, + \ 



3.3. Remarks on the curved interface case 

In this case, we consider Maxwell's system above the surface given by T : x 3 = 
(f>(x ), where x — (xi 1 X2) 1 and <j> S VF 2 (]R 2 ,R) is a scalar function. We consider 
this system in time harmonic form, in the high frequency limit, and we consider a 
perfect boundary condition on T . Again, we rewrite this system as a symetric one 

3 



iw 
e 



A°(x)u £ + Y / A^+ Cu £ - f £ (x) + A 3 u £ (x' } 0) 5. 



dx 



x 3 =(j>(x') 



(3.111) 



with C{x) given by (1.10), and f £ E L 2 (M 3 ) 3 and A" ,A> are given in (1.7), (1.8). 
We shall reduce this curved case to a plane one, by introducing the new coordinates 

y = x ,z = x 3 - (j>{x) ,x = (y,z). (3.112) 

Extending when necessary by zero in the all space M 3 , and thus (3.111) becomes 



— A°(y,z)v £ (y,z)+J2A J 

£ j=o 

where 



~ Aj dv £ {y,z) 



dv 



+Cv £ (y,z) = f £ (y,z)+A b v £ (y', 



1° 



e(y,z)Id 
r](y,z)Id 



,C = 



<r(y, z)Id 











; =o -(3.113) 



(3.114) 



are 3x3 matrix valued smoth functions, with e(y, z) , r](y, z) , <r(y, z) smooth func- 
tions in C 1 ' 



I 



I 
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Set 

v e £ (y,z) = 6(y,z)v £ (y,z) (3.115) 
and the matrix of dispersion 

3 

l(x, k) = Y,(( A )) _1 M j ( 3 - 116 ) 

where is a test function of compact support that is equal to one on a set compact 
K. 

Thus Maxwell system can be rewritten, with the cutoff function, as 

icoA°(y, z)v°{y, z) + eit ^'^f^ ? ]T A'^^-v e (x) + eC(y, z)v s £ (y, z) 

= ef B £ {y, z) + eA b v e e (y, 0) <g> 5 Z = ■ 

Then we can follow exactly the same steps as in the flat case. 
Acknowledgements: 
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